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Q^l I present a short review of the present status of the Standard Model prediction of the 
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^ ■ contributions. 
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1 Introduction 



The (yf-factor of the muon is the quantity which relates its spin s to its magnetic 
moment fl in appropriate units: 

eh 

f^ = 9ij.7, and = 2(1 + a^) . (1) 

Dirac 

In the Dirac theory of a charged spin-1/2 particle, g = 2. Quantum Electrodynamics 
(QED) predicts deviations from the Dirac prediction, because in the presence of an 
external magnetic field the muon (electron) can emit and reabsorb virtual photons. 
The correction to the Dirac prediction is called the anomalous magnetic moment. 
It is a quantity directly accessible to experiment [|. 

The experimental world average, at the time of the La Thuile meeting, 
which included the BNL published result 0) based on the 1999 yU"*" data, was 

a^(exp.) = 11 659 202.3(15.1) x 10"^° [1.3 ppm] . (2) 

There is a recent new result from the BNL collaboration II), based on /i"*" data 
collected in the year 2000, 

a/,(exp.) = 11 659 204(7)(5) x 10"^° [0.7 ppm] . (3) 

With this result, the present world average is now 

a^(exp.) = 11 659 203(8) x 10"^° [0.7 ppm] . (4) 

In this talk, I shall present a review of the various contributions to in 
the Standard Model, with special emphasis in a recent evaluation of the dominant 
contribution from the hadronic light-by-light scattering © , which has the merit to 
have stopped an avalanche of theoretical speculations, at least temporarily. 

2 Some Remarks on the QED Contributions 

In QED, the Feynman diagrams which contribute to at a given order in the 
perturbation theory expansion (powers of ^), can be classified in four classes: 

"'^See e.g. ref. for a simple and lucid review where references to the experimental literature 
can also be found. 



Diagrams with virtual photons and muon loops 

Examples of that are the lowest order contribution in Fig. 1 and the two loop 
contributions in Fig. 2. In full generality, this is the class of diagrams which 
includes those with only virtual photons, and the ones where the internal 
fermionic lines are of the same flavour as the external line. Since is a 
dimensionless quantity, these contributions are purely numerical, and they are 
the same for the three charged leptons: I = e , iJ ,t. Indeed, a] from Fig. 1 is 
the celebrated Schwinger result H) 




Fig. 1 Lowest Order QED Contribution 




Fig. 2 QED Contribution at the Two Loop Level 
while a^f^ from the seven diagrams in Fig. 2 gives the result 




At three loops there are 72-Feynman diagrams of this type which contribute. 
Quite remarkably, they are also known analytically iD. They bring in tran- 
scendental numbers like C(3)) the Riemann zeta-function of argument 3, and 
of higher complexity. 

At the four loop level, there are 891 Feynman diagrams of this type, and their 
numerical evaluation is still in progress © . 



Vacuum Polarization Diagrams from Electron Loops 

The simplest example is the Feynman diagram in Fig. 3, which gives a contri- 
bution 



6/ \2/ rrip 6b \ m,, , 
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(7) 



These contributions are enhanced by QED short-distance logarithms of the 
ratio of the muon mass to the electron mass, and are therefore very important. 
As shown in ref. , they are governed by a Callan-Symanzik type equation 



m, 
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da 
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where I3{a) is the QED-function associated with charge renormalization, and 
a^°°''(^,a) denotes the contribution to from powers of logarithms and 
constant terms. This renormalization group equation is at the origin of the 
simplicity of the result in Eq.(|^). The factor 2/3 in front of log ^ comes 
from the first term in the /3-function and the factor 1/2 is the lowest order 
coefficient of a/vr in Eq.(|^), which fixes the boundary condition to solve the 
differential equation in at the first non-trivial order in perturbation theory 
i.e., 0{^f. 




Fig. 3 Vacuum Polarization contribution from a Small Internal Mass 



Knowing the QED /5-function at three loops and (from the universal class of 
diagrams discussed above) also at three loops, allows one to sum leading, next- 
to-leading, and next-to-next-to-leading powers of log m^/me to all orders 
in perturbation theory. Of course, these logarithms can be reabsorbed in a 
running fine structure coupling a{m^). It is often forgotten that the first 
experimental evidence for the running of a coupling constant in quantum field 
theory comes in fact from the anomalous magnetic moment of the muon in 
QED, well before QCD and well before the measurement of a{Mz)- 



Vacuum Polarization Diagrams from Tau Loops 

The simplest example is the Feynman diagram in Fig. 4 below, 




Fig. 4 Vacuum Polarization contribution from a Large Internal Mass 



which gives a contribution 



aim, 



1 /rriuY hKi . rrir 
V mt m,, 



(9) 

45 Km^-J ' °'"^m/J 

In full generality, internal heavy masses in the vacuum polarization loops 
(heavy with respect to the external leptonic line) decouple. 

Light-by-Light Scattering Diagrams from Electron Loops 

It is well known that the light-by-light diagrams in QED are convergent, (once 
the full set of gauge invariant combinations is considered). Because of that, it 
came as a big surprise to find out that the set of diagrams in Fig. 5, 
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Fig. 5 Light-by-Light Scattering contribution from a Small Internal Mass 



when the lepton in the loop has a mass smaller than the external leptonic line, 
produces a contribution proportional to log(m^/me); and, in fact, with a large 
coefficient 



a 



(3)1 

ll.byl. 



2 2, rUf, 
■vr m — - 



+ ■ 



20.947. 



(10) 



This contribution is now known, analytically, for arbitrary values of the lepton 
It can be understood in the framework of effective field theories 



masses 



with different scales of masses. Again, it can be easily shown that internal 
leptonic loops with a heavy mass compared to the external leptonic line, de- 
couple. 

Altogether, the purely QED contribution to the muon anomalous magnetic moment 
of the muon, including e , fi , and r lepton loops is known to an accuracy which is 
certainly good enough for the present comparison between theory and experiment 

a^(QED) = (11 658 470.57 ± 0.29) x 10"^° . (11) 

This is the number one gets, using the determination of the fine structure constant 

a-^ = 137.035 999 59(52) [3.8 ppb] , (12) 

which follows from the comparison between the experimental determination of the 
electron (positron) anomalous magnetic moments and the QED theoretical pre- 
diction (see e.g. ref. and references therein) . The error in Eq. (|ri|) is mostly due 



to the experimental errors in the determination of the ratios of lepton masses and 
the error in the numerical integration of some of the four-loop contributions. 

The question which naturally arises is whether or not the discrepancy be- 
tween the experimental numbers in Eqs.(|^) and (^ on the one hand and the QED 
contribution from leptons alone, can be understood in terms of the extra hadronic 
and electroweak contributions predicted by the Standard Model. 

3 Hadronic Vacuum Polarization 

This is the contribution illustrated by the diagram in Fig. 6, with the shade in the 
vacuum polarization indicating hadrons. 




Hadrons 



Fig. 6 Hadronic Vacuum Polarization Contribution 



All the estimates of this contribution are based on the spectral representation 

a f°° dt 1 ^ ^ r'^ , x'^(l - x) , , 

J-M(i)X ^ ' . ■ (13) 



(h. v.p.) 



m 2 



with 

Air'^a 1 

a(t)e+e-^hadrons = — Imn(t) . (14) 

t IT 

The integration kernel in Eq.(0) shows well the underlying physical features. 

• The spectral function Imn(t), which is positive, is modulated by a known 
function of t which is also positive and monotonously decreasing. The integral 
is therefore positive and dominated by the low-energy region; mostly by the 
p resonance. 

• In QCD, the spectral function at large-t goes to a constant, which ensures the 
UV convergence of the integral. Perturbative QCD fails, however, to reproduce 
the observed hadronic shape of the spectral function below t ~ 1.5 GeV^. In 
fact, perturbative QCD with massless u and d quarks gives an IR-divergent 
result, pointing out the importance of non-perturbative effects. 



There is a lower bound to the integral in Eq.(13) 



^(..v.p.) >__^iy^ dt-lmUit), (15) 

which is governed by the slope of the hadronic vacuum polarization at the 
origin HI) . 

I have compiled in Table 1 the most recent evaluations of the integral in 
Eq.(|l3D at the time of the La Thuile meeting, and I refer to the original literature for 
the details of the various evaluations. While I was writing this talk, there appeared 
a new detailed evaluation 23), which uses the recent e^e~ data from the CMD-2 



detector at Novosibirsk, as well as the final analysis of hadronic r-decay from the 
ALEPH and CLEO detectors at LEP. Unfortunately, the results found for a'-^- " p-^ 
from the e"'"e--based data and from the r-based data are inconsistent with each 
other, even after applying radiative corrections and isospin corrections: 

(h. v.p.) _ / (684.7 ± ± 3.6,,., ) x 10-^° [e+e" - based] , 

1 (701.9 ± 4.7,.p ± 1.2,,,d ± 3.8sc,(2) ) x 10-^° [r - based] . ^ ^' 



Table 1: Compilation of recent estimates from Hadronic Vacuum Polarization 



Authors 


Contribution to x 10^° 


Davier-Hocker 
Jegerlehner 113) 
Narison ^'^'^ 
de Troconiz-Yndurain ''■^) 


692.4 ±6.2 
697.40 ± 10.45 
702.06 ±7.56 

695.2 ±6.4 



Higher order hadronic vacuum polarization contributions were first esti- 
mated in ref. . The most recent evaluation in 



gives the result 



(h.o.-h. v.p.) 



-10.0 (0.6) X 10 



-10 



(17) 



Concerning this evaluation, one should realize that there may be here a potential 
problem of double counting with some of the lowest order estimates . This is because 
part of the radiative hadronic corrections of the type indicated by the diagram in 
Fig. 7 below have already been included in the experimental cross section used to 
evaluate the lowest order contribution in Fig. 6. 




Hadrons 



Fig. 7 Higher Order Hadronic Vacuum Polarization Contribution 
This issue of double counting is under investigation at present 



H) 



4 Hadronic Light -by— Light Scattering 

These are the contributions illustrated by the diagrams in Fig. 8 below 
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Fig. 8 Hadronic Light-by-Light Contributions 



All the estimates of these contributions made so far are model dependent. There has 
been progress, however, in identifying the dominant regions of virtual momenta, and 
in using models which incorporate some of the required features of the underlying 
QCD theory. The combined frameworks of QCD in the l/A'^c-expansion and of chiral 
perturbation theory have been very useful in providing a guiding line to possible 
estimates. 



Recent progress in this domain has come from the observation EH) that, 
in large-A'^c QCD and to leading order in the chiral expansion, the dominant con- 
tribution to the muon g — 2 from hadronic light-by-light scattering comes from the 
contribution of the diagrams which are one particle (Goldstone-like) reducible; these 
are the diagrams in Fig. 9 below. The first of these diagrams (Fig. 9a) produces a 
log^ (/i/m)-term with a coefficient which is an an exact QCD result: 



log ( — ) + K{fi) 



(18) 



Here, denotes the pion coupling constant in the chiral limit {F-j^ ~ 90 MeV); the 
/i-scale in the logarithm is an arbitrary UV-scale, and m an infrared mass (either 
or m^r). 

X 7X0 





(a) 

Fig. 9 One Coldstone Reducible Diagrams in Chiral Perturbation Theory 

The dependence on the /i-scale in Eg . (p!8[) would be removed, if one knew the terms 
linear in log /j, from Fig. 9b, as well as the constant fi;(/x) from the local counterterms 
generated by Fig. 9c. Unfortunately, the determination of some of the coefficients of 
the log /i-terms cannot be made in a completely model independent wayQ neither 
the determination of the constant k{^). Nevertheless, Eq.(|T8|) plays a fundamental 
role in fixing the overall sign of the hadronic light-by-light scattering contribution to 
the muon g — 2. In the various hadronic model calculations of this contribution, there 

^Ref. H) provides a discussion of this point using a renormalization group approach. EssentiaUy 
the same arguments have been recently emphasized in ref. . 



appear indeed hadronic scales (usually the p-mass), which act as an UV-regulator, 
and play the role of fi in Eq.(|18|). Therefore, letting the hadronic scale become 
large, and provided that the model incorporates correctly the basic chiral properties 
of the underlying QCD theory, must reproduce the characteristic universal log^(/i) 
behaviour of Eq.(ff3), with the same coefficient. This test, when applied to the 



most recent existing calculations cS) (prior to the Knecht-Nyffeler calculation in 
ref. ® ) failed to reproduce the sign of the coefficient of the log^(/i)-term in Eg. (p!8|) , 
though the results from the calculations, when extrapolated to large UV-scales, 
agreed in absolute value with the coefficient of the log^(/i)-term. The authors of 
refs. 113' HD have later found mistakes in their calculations which, when corrected, 
reproduce the effective field theory test. Their results, now, agree with the Knecht- 
Nyffeler calculation @) which we report on next. 

The Knecht-Nyffeler Calculation 

In full generality, the pion pole contribution to the muon anomaly has hadronic 
structure, as represented by the shaded blobs in Fig. 10 below. The authors of 
ref. iP have shown that, for a large class of off-shell 71^77 form factors (which 
includes the laige-Nc QCD class), the contribution from these diagrams has 
an integral representation over two euclidean invariants Q\ and Q\ associated 
with the two loops in Fig. 10: 



poo roo 

/ dQl dQl W{QlQl)n{QlQl) , (19) 



where 'W{Q\,Q2) is a skeleton kernel which they calculate explicitly, and 
'H{Q\,Q\) is a convolution of two generic jF^o^,^*(fc^, /c^) form factors. In 
Large-A/^, QCD, 

J•,o^*^*(A;l^ kl)\^^^^ = (^2 _ Mf){kl- Mf) ' 
with the sum extended to an infinite set of narrow states. 

X 
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Fig. 10 Hadronic Light-by-Light from a 7r° State 



In practice, the calculation in □/ has been done by restricting the sum in 
Eq. (pOD to one and two vector states, and fixing the polynomial Cijikl^kl) 
from general short-distances and long-distances QCD properties. This way, 
they obtain the result 

^(^-oi.b.i.) ^ (5.8 ±1.0) X 10-^°, (21) 

where the error also includes an estimate of the hadronic approximation. Fur- 
ther inclusion of the rj and rjl states results in a final estimate 



(.o+r,+r,/ 1. by 1.) _ (g 3 ^ ^ 2) X 10-^° . (22) 



A Remark on the Constituent Quark Model ( CQM) 

This is perhaps a good place to comment on an argument which is often used in 
favor of the constituent quark model as a simple way to estimate the hadronic 
light-by-light scattering contribution to the muon g — 2. Since the argument 
has even appeared in print in some recent papers, I feel obliged to abrogate it 
here, so as to stop further confusion. 

The constituent quark model contribution from the diagram in Fig. 11 below. 







1^ \ < 
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Fig. 11 Hadronic Light-hy -Light in the Constituent Quark Model 



can be easily extracted from the work of Laporta and Remiddi in ref. 0) , with 
the result 



^(CQM) ^ A«\ ^ 2 

^ \nj 9 



2^^ ^ 16. 



m,. 



Mr 



+0 



m, 



Mr 



log' 



Mr 



m, 



0.616 



■ (23) 



Seen from a low energy effective field theory point of view, the constituent 
quark mass Mq in the CQM should provide the UV-regulating scale. However, 
the model is not a good effective theory of QCD and, therefore, it fails to 



reproduce the characteristic QCD log^ Mq behaviour when Mq is allowed to 
become arbitrarily large; in fact the CQM result above, decouples in the large 
Mg-limit. The argument of a positive contribution based on the CQM is 
certainly a simple argument, but unfortunately it is wrong. 

Notice however that, contrary to the naive CQM, the constituent chiral 
quark model of Georgi and Manohar 0) (see also ref. 0)) does indeed reproduce 
the correct log^ Mq behaviour in the Mq — > oo limit. This is because, in this model, 
the Goldstone particles couple with the constituent quarks in a way which respects 
chiral symmetry, and the pion pole diagram appears then explicitly. (The same 
happens in the Nambu-Jona-Lasinio model as well as in its extended version, the 
ENJL-model 0)). These models, however, suffer from other diseases^, and therefore 
are not fully reliable to compute the hadronic light-by-light scattering contribution. 
Hopefully, they will be progressively amended so as to incorporate further and fur- 
ther QCD features, in particular the short-distance constraints, following the line 
discussed in refs. ^ as already applied to the evaluation of the one particle 
(Goldstone-like) reducible diagrams in ref. ® reported above. This is why, at the 
moment, one can only claim to know the hadronic light-by-light scattering contri- 
bution with a cautious generous error, which takes into account these uncertainties. 
While awaiting for further improvement, the value quoted (at present) by our group 
in Marseille, based on the combined work of refs. 0' (appropriately corrected) 
and ref. is 

"ladlic = (8 ± 4) X 10"^° . (24) 



5 Electroweak Contributions 



The contribution to the anomalous magnetic moment of the muon from the elec- 
troweak Lagrangian of the Standard Model, at the one loop level, originates in the 
three diagrams of Fig. 12 below. 




Fig. 12 Weak Interactions at the one loop level 
^See e.g. the discussion in ref. 



where we also indicate the size of their respective contributions. Their analytic 
evaluation gives the result 0) 



,(W) 



Gf ^1 



1 



-(l-4sin^M+0f^log4 
3 V Mi mf 



m. 



Ml 



1 , 2x^(2 -X 

dx- 
1 



19.48 X 10 



-10 



(25) 



Notice that the contribution from the Higgs decouples and is very small. 

Let us recall that the present world average experimental error in the de- 
termination of the muon anomaly is 0) ^a^\^^^= ±8 x 10~^° , and, hoping for a 
continuation of the BNL experiment, it is expected to be further reduced. A theo- 
retical effort on the evaluation of the two-loop electroweak corrections is therefore 
justified. It is convenient to separate the two-loop electroweak contributions into 
two sets of Feynman graphs: those which contain closed fermion loops, which we 
denote by a^'^*^^^(ferm), and the others which we denote by a^^^^^(bos). In this 
notation, the electroweak contribution to the muon anomalous magnetic moment is 



(bos)+aJ^(')(ferm). 



(26) 



We shall review the calculation of the two-loop contributions separately. 



Bosonic Contributions 

The leading logarithmic terms of the two-loop electroweak bosonic correc- 
tions have been extracted using asymptotic expansion techniques. In fact, 
these contributions have now been evaluated analytically, in a systematic ex- 



pansion in powers of sin^ 6w, up to C[(sin^ 9 



where log 



terms, log 



terms. 



— — log T7T- terms 



Mi 



Mi 
' Ml 



terms and constant terms are kept ^) . Using 



sin^ Ow = 0.224 and Mh = 250 GeV , the authors of ref. find 



ar" '"^(bos) 



ml 



a 



v/2 87r2 71 



X 



-5.96 log 



M, 



w 



mf 



+ 0.19 



v^87r2 



a 



TT 



X 



-79.3) . 
(27) 



Fermionic Contributions 

The discussion of the two-loop electroweak fermionic corrections is more del- 
icate. Because of the cancellation between lepton loops and quark loops in 



the electroweak U{1) anomaly, in the diagrams in Fig. 13, one cannot separate 
hadronic from leptonic effects any longer. In fact, as discussed in refs. 113' 1^), 
it is this cancellation which eliminates some of the large logarithms which 
were incorrectly kept in a previous calculation in ref. It is therefore ap- 
propriate to separate the two-loop electroweak fermionic corrections into two 
classes: one is the class arising from Feynman diagrams like in Fig. 13, with 
both leptons and quarks in the VVA-triangle, including the graphs where the 
Z lines are replaced by lines, if the calculation is done in the ^^-gauge. 
We denote this class by a^^^'^\l, q) . The other class is defined by the rest of 
the diagrams, where quark loops and lepton loops can be treated separately, 
which we call a^'^^^-' (ferm-rest) i.e.. 





Fig. 13 Two Loop Electroweak Diagrams generated by the ■yjZ -Triangle 



The contribution from a^'^*^^^ (ferm-rest) brings in ml /My/ factors. It has been 
estimated, to a very good approximation, in ref. © with the result, 

Gf ml a 



a^'^^^-* (ferm-rest) 



V2 87r2 TT 



X 



-21 ± 4) 



(2^ 



Concerning the contributions to a^'^^'^\l, q), it is convenient to treat the contri- 
butions from the three generations separately. The contribution from the third 
generation can be calculated in a straightforward way, with the result 0' 



a 



EW(2). 



TT 



-3 log 



+ O 



Ml ml' 



Ml 

ml 



log 



Ml 



m7 



mt 



+ 3 



ml a 



V2 8 



TT^ TT 



X (-30.6). 



(29) 



As emphasized in ref. 113), an appropriate QCD calculation when the quark 
in the loop of Fig. 1 is a light quark should take into account the dominant 
effects of spontaneous chiral symmetry breaking. Since this involves the u and 
d quarks, as well as the second generation s quark, it is convenient to lump 
together the contributions from the first and second generation. A recent eval- 
uation of these contributions 0), which incorporates the QCD long-distance 
chiral realization as well as short-distance constraints, gives the result 



a 



EW(2). 



e, /i, u, d, s, c) 



V2 8' 



X 



-28.5 ± 1.8) 



(30) 



Putting together the numerical results in Eqs. 
Eg . (pOD , we finally obtain the value 

Gp ml 



(p8|) , (P9|) with the new result in 



V2 87c^ 



5 1 
3^3 



1-4 sin^ e 



w 



a 



TT 



(159±4) 



[15.2 ± 0.1) X 10^^°, (31) 



which shows that the two-loop correction represents a sizeable reduction of the 
one-loop result by an amount of 22% . 



6 Summary of the Standard Model Contributions 

The situation, at present, concerning the evaluation of the anomalous magnetic 
moment of the muon in the Standard Model, can be summarized as follows: 

• Leptonic QED contributions 

aQpe(/i) = 11 658 470.57 ± 0.29 x 10~^° 



• Hadronic Contributions 

— Hadronic Vacuum Polarization 

It is clear that, given the present experimental accuracy in Eq.(^), we need 
now a better understanding of the hadronic vacuum polarization contri- 
butions. Issues like the possible double counting already mentioned, and 
the improvement in the treatment of radiative corrections and isospin 
corrections 0) have now become extremely important. For refer- 
ence, I shall choose the two results from the most recent determination in 
Eq. ([TE|), combined with the higher order vacuum polarization estimate 
in Eq. (0). 



— Hadronic Light-by-Light Scattering 

(light by light) _ (8 ±4) xlO-^° 

hadronic V / 

work in progress 

• Electroweak Contributions 

Oew = (15.2 ±0.1) X 10"^° 

The sum of these contributions, adding experimental and theoretical errors in quadra- 
ture, gives then a total 

SM _/ (11 659 168.5 ±8.1) X 10-1° [e+e" - based ] , 

~ j (11 659 185.7 ± 7.4) X 10-1° [r - based ] . ^ ' 

to be compared to the experimental world average in Eq. (^ 

a;fp = (11 659 203 ± 8) x lO^^^ . 

Therefore, with the input for the Standard Model contributions discussed above, 
one finds: 



(34.5 ± 11.4) X 10-1° 3.0o- discrepancy [e+e" - based 
(17.3 ± 10.9) X 10-1° l.Qa discrepancy [r - based 



We should be prepared for a new avalanche of theoretical speculations! 
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